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Abstract 

From Noether's theorem we know symmetries lead to conservation laws. What is left to nature 
is the ordering of conserved quantities; for example, the quantum numbers of the ground state. In 
physical systems the ground state is generally associated with 'low' quantum numbers and symmet- 
ric, low-dimensional irreps, but there is no a priori reason to expect this. By constructing random 
matrices with nontrivial point-group symmetries, I find the ground state is always dominated by 
extremal low-dimensional irreps. Going further, I suggest this explains the dominance of J = 
g.s. even for random two-body interactions. 
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One of the most 'beautiful' results in mathematical physics is Noether's theorem. In 
classical mechanics a (continuous) symmetry leads to a conserved quantity, for example 
translational invariance leads to conservation of momentum, invariance under displacement 
in time leads to conservation of energy, and rotational invariance leads to conservation of 
angular momentum. 

In quantum mechanics we can think about symmetries in terms of groups and their irreps. 
Group irreps divide up a Hilbert space into subspaces; if a Hamiltonian is invariant under 
a symmetry, meaning it commutes with the generators of a group, then the Hamiltonian 
becomes block-diagonal in the irreps. 

For continuous symmetries, we label the irreps by quantum numbers, which in turn arise 
from the eigenvalues of the Casimir(s) of the symmetries. For discrete symmetries the 
Hamiltonian is still block-diagonal in the irreps, although lacking a Casimir there may not 
be a 'natural' quantum number labeling the irreps. 

While Noether's theorem is a powerful result, it tells us nothing about the relative ordering 
of states. In natural systems, however, we observe an ordering so striking and ubiquitous we 
tend to take it for granted, namely that the ground state and states lying low in the spectrum 
tend to belong to 'small' quantum numbers with the most symmetric wavefunctions. For 
example, under translational symmetry the ground state has momentum p — 0, under 
rotational symmetry the ground state has J = 0, etc. 

Of course this arises because of physics: in most Hamiltonians the kinetic energy terms 
are quadratic in linear momentum p and/or angular momentum J. This in turn occurs 
because Nature, or physicists, prefer almost-local theories, and the first nontrivial non-local 
terms are quadratic in the gradient. 

But surprisingly, even when one removes all trace of 'physics' the pattern remains. In 
nuclear structure this is seen in the discovery that many-body systems with rotationally 
invariant but otherwise random two-body Hamiltonians nonetheless tend to have ground 
states with J = 0, just like 'realistic' interactions, even though such states are a small 
fraction of the total space jl- 3]. This phenomenon is robust; for example while textbooks 
traditionally ascribe the J = g.s. to the pairing interaction {4, 5], J = still dominates 
the ground state even when the pairing matrix elements are all set to zero [6fl. It also occurs 

□ n 

in boson systems |7J such as the Interaction Boson Model[8[. Over the past decade there 
have been many proposed 'explanations.' As the distribution of many-body systems with 
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FIG. 1: Illustration of discrete rotation symmetry in the plane, C n . 
two-body interactions tend to have a Gaussian distribution of states 9| , a number of authors 



have focused on widths 10|,QJj, while others have statistically averaged in a single j-shell the 
coupling of multiple angular momenta |[12|]. As a recent Letter stated, 'the simple question 
of symmetry and chaos asks for a simple answer which is still missing [3].' 

To investigate this phenomenon, I propose a novel approach. In physics (as in art, that 
arbiter of 'beauty') we often delve deeper by stripping away assumptions to see what re- 
mains. Previous studies of random two-body interactions used shell-model diagonalization 
codes, but instead of carefully calculated matrix elements of two-body interactions they used 
random numbers, insisting only on rotational invariance. Now I go a step further: I abandon 
the shell-model framework, take random matrices, often used to investigate statistical prop- 
erties of complex systems HQ , and impose discrete rotational symmetries of the regular 
polyhedra upon them. 

Let's start with C n symmetry, the symmetry of discrete rotations in a plane. Figure [j] 
illustrates with n loci evenly spaced in a circle. The generator of discrete rotations is 

^0 ... 1 N 
1 ... 

1 ... (1) 

\0 ... 1 Oy 

which send 1 — > 2, 2 — > 3 and so on. One can easily show by hand that the most general 
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real Hamiltonian invariant under discrete rotations H = THT 1 is 
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where a,b,c, . . . are random numbers. (Note: here and throughout I do not consider adding 
reflection symmetries, saving that for future work.) 

One can solve this exactly without diagonalization by discrete Fourier transforms. One 
writes the matrix element i/jj in terms of a real function of a single index, H^j = F^j, so 
that, in ([2]), F Q = a, Fi = b,F 2 = c, etc.. The function F_j = Fj from hermiticity, and by 
further inspection one gleans F n _j = Fj + \. 

The next step is the Fourier decomposition of Ff 

i n-l 
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h k exp(i2irkj/n). 
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but using F n . 
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hk = Fj exp(— i2nkj/n) 

j=0 



[n/2] 



(3) 



(4) 



h k = F + ^2 F n cos(2nk j/n) = F + ^ F n ( j>n cos(27rkj /n) 



(5) 



where [x] is the floor function and Cj,n = (2 — 5j, n /2+i) prevents double-counting when n 
is even. It is easy to verify that h k = h n _ k . The h k are the eigenvalues of H, and one 
is justified in using k to label the different irreps. One can also find the eigenvectors ip^ k ' 
easily; k = is the most symmetric, with ip^ oc constant, while k = n/2 is arguably the 
least symmetric, with ip^ oc (— 

While © gives the eigenvalues from a simple sum, we consequently have no a priori way 
of identifying the ground state energy which irrep it belongs to. 

Next I assume there are additional degrees of freedom, as yet unspecified, and replace 
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the scalars a, b, c, . . 



bymxm real symmetric matrices A, B, C, . . . so that 

/abcd...cb\ 

B A B C ... D C 
C B A B ... E D 
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a matrix of dimension mn x mn. 

This matrix can no longer be immediately solved. What we can do, however, is to put 
H into block-diagonal form, that is, 

/ h 0...0 ^ 
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where 
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h fe = F + F i Cj>cos(27r/cj/n) 



(8) 



n=l 



Here comes the key step. While we cannot analytically compute the ground state energy 
of each h^, we can compute the variance. If we assume the Fj are independent random 
matrices, each with the same variance <Tq (assumptions which can be relaxed), then the 
variance of the kth block is 



[n/2] + l 

** = *o(l+ E (l n co S 2 (27rk(l-j)/n) 
i=2 



(9) 



For k = 0, this yields approximately a^(2n-\-l), while otherwise this will yield approximately 
(Xpfn + 1) (because the average of cos 2 x ~ 1/2). Thus the matrices with k = will have the 
larger widths and the ground state k will be one of those two. 

Now C n is an abelian group, but we can do the same analysis for other, nonabelian point- 
symmetry groups based on regular polyhedra. First, consider the tetrahedron, illustrated in 
Fig. [2j The most general Hamiltonian invariant under any discrete rotation about any of its 
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FIG. 2: Illustration of tetrahedron space symmetry. 



facets is 
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Bringing this into block diagonal form, 
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(10) 



A + 3B 
A- B 
A- B 
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(11) 



That is, there is an irrep of dimension 1 with matrix A + 3B and an irrep of dimension 
3 with matrix A — B. Again assuming A, B are independent but have the same variance 
the variance of the 1-dimensional irrep (which has the most symmetric eigenvector) is IOctq 
while that of the 3-dimensional irrep is 2<7q. All else being equal, the ground state is much 
more likely belong to the 1-dimensional irrep. 

With the basic idea in hand, it is easy to consider other polyhedra. For the octahedron, 
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FIG. 3: Illustration of octahedral space symmetry. 



Fig. [31 the invariant hamiltonian is 

^ABCBBB^ 
B A B C B B 
C B A B B B 
B C B A B B 
B B B B A C 
B B B B C A 

which has a 1-dimensional irrep, with matrix A + 4B + C, a 2-dimensional irrep with matrix 



V 
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A — 2B + C, and a 3-dimensional irrep with matrix A — C. The variances are 18a^, 6<Jq, 
and 2<7q, respectively. 

For the cube, Fig. HJ the Hamiltonian is of the form 

^ A B C B B C D 
BABCCBCD 
CBABDCBC 
BCBACDCB 
BCDCABCB 
CBCDBBBC 
DCBCCBAB 
CDCBBCBA 
with two 1-dimensional irreps with matrices A ± 3B + 3C ± D and two 3-dimensional irreps 
with matrices A — C ± (B — D). The 1-dimensional irreps have variance 20<Jq while the 
3-dimensional irreps have variance 4ctq. 



V 
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FIG. 4: Illustration of cubic space symmetry. 



While it remains to be proved in general, the lesson is clear: Starting with random 
matrices and imposing symmetries, the ground state is naturally dominated by certain irreps, 
generally irreps with lowest dimension. 

Now I turn to continuous symmetries such as SU(2). I consider wavefunctions that can 
be written in the form ipi(v)Yi m (Q , 0) where all of the rotational information is bound up in 
the spherical harmonic Yi m and v refers to internal degrees of freedom [15(. I'll return to 
the latter in a moment. Using the angles from spherical coordinates, the Hamiltonian is of 
the form H(9'<j>',9</>), but imposing rotational invariance means H can only depend on the 
angle u between 6', <$' and 9, <p as given by cos u = cos 9 cos 9' + sin 9 sin 9' cos(0 — <p'). Then 
H{9'(j)' , 9(f) = F{oj), where F(u) = F(u> + 2tt) is a periodic function and, using Hermiticity 
(and assuming H is real) F(u) = F(2ti — uj) is symmetric with respect to u = tt. Expanding 

2 7+1 

F(cu) = J2hj^ r Pj(cosu) 
= H h J E YjuWiMYjM&fi (14) 

J M=-l 

so clearly the hj are again the eigenvalues, with Yjm as eigenfunctions and with the eigen- 
values independent of N, as one expects. 

Once more I assume F(u>) to be a matrix-valued function, and thus hj to be a symmetric 
matrix given by 

P'K 

hj = 2tc / Pj(cosu)F(u) smudu. (15) 
Jo 

As before, let a be the variance of the matrix elements F independent of u. Then I estimate 
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the variance of the matrix elements of hj 

a 2 = An 2 a 2 jT Pj (cos u) sin 2 udu (16) 

Eqn. (fl6|) can be computed numerically, and leaving off the factor An 2 a 2 , the values are 
1.571, 0.393, 0.245, 0.178, 0.139 for J = 0, 1, 2, 3, 4 respectively. 

This suggests that in many-body system, subspaces with low-valued quantum numbers 
will have larger widths. But in realistic, finite many-body calculations, subspaces with 
different Js have different dimensions. Furthermore, in the above argument each hj has 
independent random matrix elements, which typically has a semi-circular density of states 



14j , yet for many-body systems with only two-body interactions the density of states tends 



towards a Gaussian [ 

I can approximately correct both deficiencies. First, following standard results on matri- 
ces with Gaussian-distributed matrix elements [ijj], I let a'f = y/T^jaj be the width of the 
subspace of states with angular momentum J. Then, for each J, I simply create Nj energies 
via a random Gaussian distribution of width a} 5 , and ultimately determine the fraction of 
ground states with angular momentum J. 

Finally, I compare against a simulation via configuration-interaction calculations, that is, 
diagonalizing the Hamiltonian for fixed numbers of particles in finite single-particle spaces. 
Figure [5] shows the case of eight fermions (neutrons) in the IP1/2-IP3/2-O/5/2-O/7/2 shell-model 
space. 

I take an ensemble of rotationally invariant, two-body but otherwise random interactions, 
and tabulate the fraction fci of states that have a given angular momentum J. This should 
be compared with the native fraction of states with that J in each many-body space, /space = 
Nj/N tot (N tot is the total dimension of the many-body space), and the fraction with J = 
is dramatically enhanced. I also compare with the fraction of ground states with a given J 
predicted by my simple random matrix picture, /rm- The only input are the dimensions Nj 
and the universal variances computed in f fl6|) . 

For such a simple picture, the random matrix model yields qualitatively excellent results, 
generally predicting the enhancement or suppression of different Js in the CI simulations 
relative to the native fractions /space- The RM model successfully predicts an enormous 
enhancement of J = in the ground state, for this case and many others not shown due to 
lack of space. 
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FIG. 5: (color online) Distribution of ground state quantum numbers for eight neutrons in the 
lp-Of shell. Empty bars are the fraction of states in the model space with a given J, horizontally 
striped (red) bars are the fraction of ground states with a given J predicted by a random matrix 
(RM) model, and vertically striped (blue) bars are fraction of ground states with a given J from 
configuration-interaction (CI) diagonalizations of an ensemble of random two-body Hamiltonian in 
a shell-model basis. 

This analysis suggests the predominance of angular-momentum zero ground states is pri- 
marily a function of the width of the angular-momentum-projected many-body Hamiltonian; 
furthermore, the width is largely decoupled from the microphysics, instead depending only 
on the projection integrand (fT6|) and on the dimensionality of subspaces with good quantum 
numbers. The simplicity and decoupling from the microphysics may be why the phenomenon 
is so robust and so universal. 

The U.S. Department of Energy supported this investigation through grant DE-FG02- 
96ER40985. 
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